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STEADY  AND  UNSTEADY  LOADINGS  AND 
HYDRODYNAMIC  FORCES  ON  COUNTERROTATING 
PROPELLERS 


S TSAKONAS.  W JACOBS  and  M AL1 

Stevens  Institute  of  Technology.  Hoboken.  New  Jersey.  USA 


The  Linearized  lifting  surface  theory  his  been  applied  to  the  evaluation  of  blade  loadings  and  re- 
sultant hjdrodyn imic  forces  and  moments  ( thrust , torque,  bearing  forces  anl  bending  moments)  of  counter- 
ntating  propeller  systems  with  equal  anl  unequal  number  of  blades  operating  in  uniform  and  nonuniform 
inflou  fields,  both  units  rotating  with  the  sane  RPM.  The  mathematical  model  takes  into  account  as  real- 
istically as  possible  the  -geometry  of  the  propulsive  device,  the  mutual  interaction  of  both  units  ani  the 
three-dimensional  spatially  varying  inflow  field.  A computer  program  Ivin  been  developed  adaptable  to  high- 
speed digital  computers  (CDC  6600-  7600 ) for  counterrotating  propulsion  systems  of  equal  and  unequal  blade 
number  in  uniform  inflow  for  comparison  with  experimental  results  if  available. 


I.  INTRODUCTION 

The  combination  of  two  counterrotating  pro- 
pellers on  fast  ships  has  been  shown  to  offer 
considerable  improvement  in  propuls i ve  efficiency 
when  compared  with  a single  screw  LlJ.  Further- 
more, since  the  total  reguired  power  is  divided 
between  two  propellers,  this  results  in  a reduc- 
tion of  blade  loadings  and  hence  the  inception  of 
cavitation  is  delayed.  These  are  the  main  advan- 
tages of  this  propulsive  system;  its  principal 
disadvantage  lies  in  the  mechanical  complications 
in  transmitting  power  through  a coaxial  counter- 
rotating shaft. 

The  CR  (counterrotating)  propulsive  system  is 
also  expected  to  have  more  favorable  vibrational 
behavior.  From  tests  of  a 4-0-5  CR  system  (4- 
bladed  forward  and  5-bladed  after  propeller)  in  the 
wake  of  a model  of  a fast  cargo  liner  l2J,  it  ap- 
pears that  the  ratios  of  amplitudes  of  excitation 
to  mean  thrust  are  comparable  to  those  of  a single 
screw  providing  the  same  power.  However  in  these 
tests  the  nonuniform  wake  is  by  far  the  dominating 
cause  of  vibration.  The  effects  of  the  interaction 
of  both  propellers  are  small  in  comparison  and  the 
higher  frequency  excitations  cannot  be  determined 
at  al I accurately. 

A better  understanding  of  the  mechanism  of 
the  interaction  can  be  obtained  by  considering  the 
CR  system  under  open-water  conditions  (uniform  in- 
flow field)  so  that  wake  harmonics  are  not  present 
to  mask  the  interaction  phenomenon. 

The  calculation  procedure  is  based  on  the 
analysis  of  Reference  3 for  the  cases  of  CR  sys- 
tems of  equal  and  unequal  blade  number,  operating 
at  equal  or  unequal  RPM  in  uniform  and  nonuniform 
inflow  fields.  In  that  reference  the  true  geom- 
etry of  the  helicoidal  blades  was  taken  into  ac- 
count with  the  exception  that  blade  thickness  was 
assumed  negligible. 

In  the  present  paper,  CR  systems  of  equal  and 
unequal  number  of  blades  will  be  considered  operat- 
ing at  equal  RPM  in  a uniform  inflow  field.  The 
blade  thickness  effects  will  also  be  considered,  as 
additional  velocity  perturbations  on  the  LH  (left- 
hand)  sides  of  the  two  surface  integral  equations 
which  state  the  kinematic  conditions  on  both  units 


of  the  CR  system. 

The  development  of  this  pair  of  surface  in- 
tegral equations  is  based  on  a linearized  unsteady 
lifting  surface  theory  as  adapted  to  the  marine 
propeller  case.  Their  kernel  functions  are  derived 
by  means  of  the  acceleration  potential  method  and 
the  surface  integrals  are  reduced  to  line  integrals 
by  employing  the  mode  approach  in  conjunction  with 
the  "generalized  lift  operator"  technique  L4J.  Then 
by  the  collocation  method  the  line  integral  equa- 
tions are  reduced  to  two  simultaneous  sets  of  al- 
gebraic equations.  Finally,  the  solution  of  these 
is  obtained  by  an  iterative  procedure,  assuming  at 
first  that  the  effect  of  the  after  propeller  on  the 
forward  propeller,  except  for  the  velocity  field 
due  to  the  thickness  of  its  blades,  may  be  neg- 
lected. 

The  computation  procedure  adapted  to  a high- 
speed digital  computer  (C0C-6600  or  7600)  will  be 
utilized  if  experimental  information  will  be  avail- 
able. 

NOMENCLATURE 

A subscript  index  of  after  propeller 


F subscript  index  of  forward  propeller 

Fx  y z forces  in  axial,  horizontal  and  vertical 
directions 

l^(x)  defined  in  Equation  (7a) 

I (v)  modified  Bessel  function  of  order  m 

i subscript  index  of  control  point 

j subscript  index  of  loading  point 

K (v)  modified  Bessel  function  of  order  m 
w 

Kj j kernel  function  of  integral  equation 

K . j , K • j modified  kernels,  after  chordwise- 
J J integrations 

L(r)  spanwise  loading  distribution,  Ib/ft 

L^(p)  spanwise  loadinq  components  (coefficients 

of  chordwise  distribution),  Ib/ft 
£.  i nteger  multiple 
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m order  of  lift  operator 

m,  index  of  summation 

k 

N^.N^  number  of  blades  of  forward  and  after  pro- 
pel lers 

n order  of  chordwise  mode 

n blade  index 

P perturbation  pressure 

Q moments  about  x-,  y-  and  z-axes 

* * y » z 

q.  order  of  blade  harmonic 

r radial  coordinate  of  control  point 

r superscript  index  of  control  point 

r^Q  after  propeller  radius 

rFo  forward  propeller  radius 

Sj  lifting  surface 

t time 

U uni  form  veloci ty 

u variable  of  integration 

V.(r)  Fourier  coefficients  of  velocity  normal  to 
the  blade 

Wjfx.r.'Pjt)  induced  velocity  at  control  point 

x.x1  longitudinal  coordinate  of  control  point 

x(x'),r,cf>  cylindrical  coordinate  system  of  control 
points 

y horizontal  Cartesian  coordinate 

z vertical  Cartesian  coordinate 

B hydrodynamic  pitch  angle 

e distance  between  the  two  propeller  planes 

d(n)  chordwise  mode 

Sj  angular  coordinate  of  loading  point 

angular  chordwise  location  of  loading  point 

^bF’^bA  Pr°jected  semichord  length,  in  radians,  of 
forward  and  after  propellers 

6.  -j^-fn-l)  • n='  • • • • 

J n Nj  J 

6 ( r)  geometric  pitch  angle 

(x)  defined  in  Equation  (7b) 
positive  integer  multiple 
5.J'  longitudinal  coordinate  of  loading  point 

5(D.P.6  cylindrical  coordinate  system  of  loading 
points 

0 radial  coordinate  of  loading  point 

0 superscript  index  of  loading  point 

Pp  mass  density  of  fluid 

ct  angular  measure  of  skewness 

4(m)  generalized  lift  operator 

tp.  angular  coordinate  of  control  point 

<pa‘  angular  chordwise  location  of  control  point 

»1  angular  velocity  of  propeller  (absolute 

value) 


II.  LINEARIZED  UNSTEADY  LIFTING  SURFACE  THEORY 

Two  counterrotating  propellers  are  operating 
in  the  flow  of  an  ideal  incompressible  fluid.  The 
propeller  arrangement  and  the  coordinate  system  are 
shown  in  F i gure  I . 

The  basic  relation  of  the  interaction  phenom- 
enon is  that  the  negative  velocities  induced  by  the 
propulsion  system  on  each  propeller  lifting  surface 
should  be  balanced  by  the  downwash  velocity  distri- 
bution at  that  surface,  thus  expressing  the  require- 
ment of  an  impermeable  boundary.  The  kinematic 
boundary  conditions  on  both  lifting  surfaces  are 
expressed  as  two  simultaneous  surface  integral  e- 
quat ions  : 

Vvw^  “JT  APF(SF.PF.6F;t) 

SF 

KFF(  W WPF'V^dSF 


+ JT%(5;,PA.eA;t) 


KAF(VVV?A'PA-eA:t>dSA 


WA(xA'rA'<pA;t)  = If  APF(5F.PF.eF;t) 
SF 


KFA(xAirA’V?F*VeF:t)dSF 


+ if  6pA<?;-pA-eA:t> 

SA 


xtx'J.r.tp  and  5(5'), P, 9;  cylindrical  coor- 
dinates of  control  and  loading 
points,  respectively 

F and  A;  subscripts  indicating  forward  and 
after  propel ler 

t:  time,  sec 

$F,S  • forward  and  after  propeller  sur- 
faces. ft2 

WF,WA : velocity  distributions  normal  to 

forward  and  after  propel lers , ft/sec 

APp.PP^:  unknown  loadings;  pressure  jumps 

across  the  lifting  surf aces . Ib/ft2 , 
i.e.,  AP=P. -P+  pressure  difference 
between  back  (suction  side)  and 
face  (pressure  side) 

K..:  kernel  function  representing  the  in* 
■*  ' duced  velocity  on  an  element  i of 
a blade  due  to  unit  amplitude  load, 
located  at  each  and  every  element  j , 
fb/lb-sec 

The  second  term  on  the  RH  (right-hand  side) 
of  Eq.(l)  and  the  first  term  on  the  RH  of  Eq.(2) 
are  the  interaction  effects.  The  remaining  terms 
are  the  self-induced  velocities  by  the  individual 
propel lers . 

The  unknown  loadings  and  the  onset  velocity 
distributions  are  cyclic  in  nature.  Then  for  a 
CR  system  with  right-hand  aft  propeller  and  left- 
hand  forward  propeller  rotating  at  equal  RPM 
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( O -i  Vrt 

AP  (Sp.Pp.b F;t)-Re  L,  aPF  (5  p 0 )e 
Xk=0 


® (O  >Mit 

<^VVe 

k (3) 


",  _ *iqcflt 

WF(xF,rF,tpF;t)»Re  L,  (xF,  rF,<fy)e  p 

qF=° 


® _ (q«)  iq.flt 

WA(XA-rA’;pA:t)cRe  L.  WA  <*;-rA'Ve 

qA=° 
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where  qj  designates  the  order  of  shaft  frequency 


or  order  of  harmonic  of  the  inflow  field. 


that 


of  the  loading  distribution  to  be  determined  by 
the  analysis,  and  ft  is  the  absolute  value  of  the 
angular  velocity  of  each  propel ler, (q;  and  ^k  are 
both  positive  integers.)  The  known  downwash  veloc- 
ities and  the  unknown  loadings  are  expressed  in 
complex  conjugate  form  in  (3)  and  (4),  where  finally 
the  real  part  is  taken. 

The  velocities  W;  are  caused  by  flow  dis- 
turbances such  as  those  due  I)  to  wake,  2)  to  in- 
cident flow  angle  which  is  the  difference  between 


blade  and  the  hydrodynamic  pitch  angle  P=tan_lll/Wr 
where  U is  forward  speed  and  r is  the  radial 
location  of  the  corresponding  helix,  3)  to  blade 
camber,  4)  to  "non-p I anar"  blade  thickness,  and  5) 
to  the  effects  of  the  thickness  of  the  blades  of 
each  propeller  on  the  velocity  field  of  the  other. 
Within  the  limits  of  the  linear  theory,  the  effects 
of  the  flow  disturbances  can  be  obtained  separately 
and  then  simply  added  together. 

Although  the  analysis  applies  to  both  non- 
uniform  and  uniform  inflow,  as  mentioned  earlier, 
the  solution  by  an  iterative  procedure  will  be  re- 
stricted to  the  uniform  inflow  case  (no  wake). 
Furthermore,  the  disturbance  due  to  the  so-called 
"non-planar"  thickness  (since  a propeller  blade 
lies  on  a helicoidal  surface  of  variable  pitch,  its 
thickness  affects  its  own  velocity  field)  will  be 
ignored  as  negligibly  small  L 5 J - 

After  the  chordwise  integrations  are  per- 
formed by  means  of  the  mode  approach  and  the 
generalized  lift  operator  technique,  the  pair  of 
surface  integral  equations  (I)  and  (2)  are  reduced 
to  the  following  set  of  line  integral  equations 
for  given  q;  , order  of  shaft  frequency,  given  m , 
order  of  lift  operator,  and  ri  , order  of  chordwise 
mode  shapes,  for  the  case  of  equal  RPM: 


wF(VS) 


(qc.fi)  ® ,-  -, 

(r  )-Jlf  (Bp)  L C ’ (m,»q  +ijN  ) dp 

Pp  171}= 


. ® “ (ka=q  -2X2N«,n) 

+J  ^ ^ LA 

Pp  *2=0  ma=C 
j(m.n) 


(PA) 


AF 


(ma=qF-*iNA)dpA 
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TT  (rA)’-’  L --F 
Op  X3”0  m3»0 


® (A.q“q  +2i3Nr,  n) 
VL,  <PF) 


, “(m.n) 


FA  (w»“Vi3MF)dpF 


VA 


A <pA,mjLRAA,B)  K"V*«NA)dPA  (6> 


VIII  JJ 


(n)  , 


Here  L' ” 7 ( Pp)  and  (pA)  are  the  unknown 

normal  loading  components  of  the  n chordwise  mode 
for  each  blade  in  Ib/ft  of  span,  Np  and  NA  are  the 
blade  number  of  forward  and  after  propellers,  and 
ik  is  integer.  The  bars  and  superscripts  m and 
n indicate  that  the  quantities  have  been  inte- 
grated along  the  chord. 

The  values  of  kk  and  ttik  shown  in  (5)  and 
(6)  are  arrived  at  by  equating  the  time-dependence 
on  LH  and  RH  sides  so  that 


-iqFnt 

-U.ftt 

e *» 

e 

for  the  first  term  of  the 
first  integral  equation 

-iqpUt 

i (\2-2m2)Ot 

e ■ 

e 

for  the  second  term  of  that 
equat ion 

i qA«t 

-i  (Xj-2m^)  fit 
e 

e » 

for  the  first  term  of  the 
second  integral  equation 

iqAnt 

i\4nt 

e 

e 3 

for  the  second  term  of  that 
equat ion 

and  from 

the  summation 

1 over  all  blades  of  a pro- 

peller  which  is  represented  by 

N ±'(mk"^k^n  rN 

for  mk-Xk=\N,ik“0,±l.  ... 

Lj  e 
n=  1 

“ lo 

otherwise 

where  6n  = 2n(n-l)/N 


The  respective  kernels  are  derived  in  Refer- 
ence 3 for  RH  forward  propeller  and  LH  after  pro- 
peller. For  IH  forward  and  RH  after  propeller, 
they  are  given  in  the  following  section  in  final  form. 


III.  THE  KERNEL  FUNCTIONS 


The  Kernels  KFF  and  i<AA 


These  functions  describe  the  self-induced 
velocity  at  a point  on  a propeller  blade  due  to 
unit  amplitude  load  at  various  locations  on  all 
the  blades  of  the  same  propeller.  The  development 
for  a right-hand  propeller  is  given  in  Reference  6 
and  yields 


£ K^m'n\m=q+iN)  = {- — — -} 

m=-»  ^4ttp,lt  r_ 


- i qAJ 


"OfWV  aZiTTT* 


k.  {.(•>- if  *44?^*} 


S=q?*N 


(7) 


where 


9(u)  = (IK)mB'(u)e 


i — A a 
a 


("On 


1 


rlm(  lu+aiNI  p)Km(  I u+aiNI  r)  for  p<  r 


l-lm(  I u+aiNI  r)Km(  I u+a£NI  p)  for  r < p 
B'(u)  = ^au  + a3  £n  + 

vu+aa*N+  X>(R\<OebpJ 

PF  = fluid  mass  density,  slugs/ft3 
r » propeller  radius,  ft 


r p 

Aa=o-o  = difference  between  skewness  of  the 


blade  at  the  control  point  r and 
skewness  at  a loading  point  p,  radians 


a “ ft  r0/U  and  p and  r are  also  non- 
d imens ional i zed  by  rQ 
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0^,0^  ■ subtended  angle  of  projected  semi  chord 
of  blade  at  r,  at  c,  radians 

lm(  ) = modified  Bessel  function  of  first  kind 

Km(  ) ■ modified  Bessel  function  of  second  kind 

l = 0,  ±1 , ±2,  ... 

In  this  equation  the  chordwise  integration  is  repre- 
sented by 

|(S)(X)  =i  J *(m)eixcos^dg>0  (7a) 


where  4(r?i)  is  the  lift  operator  function  and 

A^(y)  - ^ J Q(n)e  'YC°Se“sin0odea  (7b) 
o 

where  c(n)  is  the  chordwise  mode  shape  selected. 
(See  References  3, A, 5, 6.) 

In  Eq.(6)  the  kernel  function 

Li  k!*’^(hf4  = q.+24N.) 
m4=-«  AA  * MA  * A' 

is  given  by  Eq . ( 7)  with  q=qft.  N=N^,  r=r^  and  P=C^. 
However,  r0~rf.Q  the  radius  of  the  forward  propeller. 
In  Eq.(5)  the  kernel  function 

mjLRirn)<mi=-qF+MF) 

is  also  given  by  Eq.(7)  but  with  q=qF,  N=Np,  r=rp 
and  e=Pp.  In  both  Eqs.(5)  and  (6),  the  radial 
positions  r and  p and  the  inverse  advance 
ratio  a are  non-dimensional i zed  by  forward  pro- 
peller radius  rpo 

The  Kernels  i<AF  and  K 

These  are  the  kernels  of  the  cross-coupling 
terms  of  Eqs.(5)  and  (6).  Let  the  distance  be- 
tween the  propeller  planes  of  the  two  units  of  the 
CR  system  be  e (in  terms  of  rpg)  . Then  for  a RH 
after  propeller  at  a distance  e from  a LH  for- 
ward propeller  operating  at  the  same  RPM,  the  de- 
rivations of  Reference  3 can  be  reduced  to  the 
following  final  forms; 

_(m,n)  ' ’na  ) rF  im  (2a  +ae) 

K'  '(ma-q  + A.N  .20)=;- V-f— L e 

AF  FA  '•bnp.ii^r  J „ T i 4 

^fl  TO  a v |+a2r2 

iq_(o_-o  -ae)  , . ® 

• e {a(0)B(0)-£  J[a(u)B(u>AHB(u)]~} 

o 


r -kna(2o  -ae) 


™,n\m3«qA+i3NF«:0)-'[— - — t •]> e 

k A F ^PfU  rFOjavro 


A(u)=| 


lm2  (lu-aftTte-q^lP^I^J  lu-a(m2-qF)|rF)  for  p^sr^. 
lnfe(lu-a(m2-qF)lrF)Km2(  lu-a(nz-qF)|PA)  for  rp<PA 


r 2 rn2n  r 2 ^2  n 

B(u)-  |^au-a  (m2-qF)  - [au-a  (ma-qf.)+  -j-J 


iu(a  -o.-ae)/a 

• F A '(m)(<-‘v->bF; 


A(n\(2ma  - qp  - 


-iqA(aF-°A-ae)f 


{c(0)D(0)  - i ;[C(u)0(u)-C(-u)D(-u)]^}  (9) 


lm3(lu+a(m3 -qA)lcF)Km3(lu+a(m3-qA)lrA)  for  PFsrft 

C(^{ 

lm3(lu+a(m3-q^lr^Km3(lu+a(m3-qA)lpF)  for  r A<pp 

D(u)=  j^au  + a2  (m3-qft)  + -^]  [au  + a2(m3-qft)  - ^ J 

-iu(a  -a  -ae)/a  ... 

•*  F * 

• A<n\(-2na  * qt  - Ji  \rJ 

The  kernels  have  been  programmed  with  proper 
consideration  being  given  to  evaluating  the  finite 
contributions  of  the  Cauchy-type  singularities  in 
the  u- i ntegrat i ons  at  u=0  and  of  the  Hadamard- 
type  higher  order  singularities  in  the  p-integra- 
t ions  when  p-»  r . 

IV.  THE  NORMAL  VELOCITIES 

The  LH  sides  of  the  integral  equations  repre- 
sent the  normal  components  of  the  velocity  pertur- 
bations above  that  producing  zero  loading  (lift) 
which  corresponds  to  a rotating  thin  plate  lying  on 
the  helicoidal  surface  of  pitch  angle  (i.e.,  hydro- 
dynamic  pitch  angie) 

-I  U 

b = tan  — 

Ur 

■/here  U = forward  speed,  r = radial  location  of 
the  corresponding  helix.  The  perturbations  con- 
sidered are  those  due  to  hull  wake  (non-uniform  in- 
flow field),  blade  camber,  incident  flow  angle,  and 
also  due  to  the  thickness  of  the  blades  which 
mainly  affect  the  velocity  field  around  both  pro- 
pellers of  the  CR  system.  The  effect  of  "nonplanar" 
thickness  which  contributes  to  the  loading  of  both 
units  is  neglected  as  being  small  L5j.  The  effects 
of  each  of  these  imposed  flows  on  the  blade  are 
calculated  separately  and  simply  added  together  as 
allowed  by  the  linearity  of  the  theory. 

The  left-hand  sides  of  the  integral  equations 
due  to  the  wake  contribution  can  be  harmonically 
analyzed  and  written  in  the  form 

W(r,3)  = £ V n(r)e  'qiJ 
q=0  n 

where  cl  can  be  expressed  in  terms  of  the  moving 
coordinator  system  attached  to  each  propeller  by 

SF  = and  yA  = “e0A + ^Ar  for  the  case 

the  forward  and  after  propeller,  respectively,  as 
shown  in  Figure  I.  The  normal  wake  velocities 
Vn(r)  can  be  determined  from  the  harmonic  analysis 
of  the  wake  measurements,  as  shown  in  Reference  5. 
After  the  trigonometric  transformation 

0o  = a - ebcosea 

and  application  of  the  "lift  operator"  of  order  m 

n 

j*  «(S)W(r,0)d9a 
o 
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the  following  expressions  result  for  the  left-hand 
sides  of  the  pair  of  integral  equations  relating 
the  unknown  loadings  with  the  qiven  "downwash"  at 
each  prope tier; 

-(qA-S)  (V 

v '>i . v»- v -;va  («),  . 

U IJ  lqA  bA' 

for  the  after  propeller,  and  (10) 


(qr.m) 

Wp  (rp) 


(qF) 

VF  (rF)  'iqF°F, 


,(qFv] 


for  the  forward  propeller.  It  should  be  noted 
that  the  factor  expiq^t  has  been  omitted  from 
the  above  expressions. 

The  velocities  induced  by  the  incident  flow 
angle  and  camber  effects  are  independent  of  time 
because  the  blades  are  considered  rigid  so  that 
only  the  steady-state  loadinqs  will  be  affected. 

For  the  flow  angle  (f)  effects,  the  dimen- 
sionless perturbation  velocities  after  the  appli- 
cation of  the  lift  operator  become: 


been  shown  to  be  a good  approximation  for  deter- 
mining velocity  and  pressure  t8,9j  on  a point  in 
the  neighborhood  of  an  operating  propeller  as  long 
as  it  is  not  a point  on  its  blade  and  particularly 
near  the  leading  edge.  It  should  also  be  recalled 
that  the  velocity  and  pressure  fields  generated  by 
an  operating  propeller  even  in  a uniform  inflow 
yields  steady  and  unsteady  components  of  the  re- 
spective field.  Thus  although  it  is  independent 
of  time,  the  blade  thickness  produces  both  steady 
and  unsteady  components  of  the  velocity  field. 

Following  the  same  procedure  as  in  Reference 
8,  it  can  be  shown  that  the  dimensionless  velocity 
normal  to  the  blades  of  the  forward  propeller  in- 
duced by  the  after  propeller  thickness,  with  maxi- 
mum thickness-chord  ratio  tQ/c  , is  given  for 


(0,m) 

(0F(rF>)  4a‘rFNA 

^ U '‘A  WiwV 


T^A> 


iu(a  -a  +ae)/a 


r v A F 

ou(IK^F(u.PA)R.P.[e 


'<7  VJdudpA  (,3> 


L%F(rF)-SF(rF)jl 


lo<upA>Ko(urF)  f°rBA<rF 
° L|0(urF)Ko(uPA)  forrF<PA 


V—J  LepA(rA>-CA(rA>],(m)(0> 


F(u,pa)  = 


UC>.  A UtjLA  * Ufa,  A . 

sin  M ^ bA  . bA 

a a \ a / 


for  the  forward  and  after  propellers,  respectively 
w here  is  the  blade  pitch  angle  and  3 is  the 

hydrodynamic  pitch  angle  of  the  reference  helicoi- 
da!  surface  (i.e.,  of  zero  loading). 

For  the  camber  (c)  effect  of  the  forward 
propeller,  the  dimens  ionless  velocity  ratio  after 
the  application  of  the  lift  operator  becomes 


(0,m) 

WF  (rp> 


df(rF.sF) 


F(rF,sF)=  camberline  ordinates  as  fraction 

of  expanded  chord  length,  measured 
from  the  face  pitch  line 

Sp  = ( I -cos'-Pq,) /2  , chordwise  location 
non-dimensional i zed  on  the  basis 
of  chord  length 

For  the  after  propeller,  the  same  expression 
is  valid  and  only  the  subscript  F must  be  re- 
placed by  A . 

The  evaluation  of  the  integral  of  (12)  is 
given  in  Reference  7 for  arbitrary  camber  shape. 


Iq(  ) and  Kq(  ) are  modified  Bessel  functions 
of  first  and  second  kind  of 
order  zero 

and  for  qF=2iNA  where  j£=I,2  ...  , by 

(2iN  ,m) 

lQ F (rF)J  „ . ',a"rFNA  -i%(2Va£) 

U tA  " v'|+aferF 

(p  ) 

P.  bA  MCA 

A 

CO 

‘ jF(u,PA)LG(u)-G(-u)]dudpA  (14) 

where 

G(u)^au+ANA(a3-y]|£  (,Ufa£NA|PA)K£  (|u^NA|r  ) 

rF  A A 

i U (O  — ° +3€)  / 3 . \ 

F .“Vv 


V.  THE  EFFECTS  OF  BLADE  THICKNESS  OF  EACH  PRO- 
PELLER ON  THE  VELOCITY  FIELD  OF  THE  OTHER 

In  addition  to  the  disturbances  of  the 
velocity  field  about  each  propeller  due  to  wake, 
its  flow  angle  and  camber,  the  disturbances  con- 
sidered are  those  due  to  the  effects  of  the  thick- 
ness of  the  blades  of  one  propeller  on  the  velocity 
’incident  on  the  other.  These  normal  velocities  on 
the  LH  side  of  the  integral  equations  have  been 
developed  on  the  basis  of  "thin"  body  approxima- 
tions L 8 ] . Furthermore,  it  is  assumed  that  the 
thickness  distribution  is  approximated  by  a 
lenticular  cross-section,  an  assumption  which  has 


f°r  P^<rF.  If  PA>  rF,  these  factors  are  inter- 
changed in  the  modified  Bessel  functions 

The  velocity  normal  to  the  blades  of  the 
after  propeller  induced  by  the  forward  propeller 
thickness  is  for  qA=0 
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where 


(15) 


(IK), 


and  for  = 2XNp,  4=1,2,  ... 
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whe  re 


»(„)-(>», l»Fl>,"-i,;]l£|^(|u,aI»Flpr)K1^(lu.alMF|rA) 


for  0F<  rft.  If  pF  > rft, 


these  factors  are  inter- 
changed in  the  modified  Bessel  functions. 


VI.  SOLUTION  OF  THE  PAIR  OF  INTEGRAL  EQUATIONS 

It  is  seen  from  Eqs.(5)  and  (6)  that  the 
loading  on  each  propeller  of  the  CR  system  is  af- 
fected by  all  the  harmonics  of  the  inflow  field. 
Here,  however,  the  solution  will  be  limited  to  the 
uniform  inflow  case,  i.e.,  to  the  steady-state 
normal  velocities  due  to  camber  (c)  and  to  inci- 
dent flow  angle  (f),  and  to  the  steady  and  unsteady 
effects  of  the  interaction  (i.e.,  cross-coupling 
terms)  between  the  propellers  and  to  their  respec- 
tive thicknesses. 

Equation  (5)  becomes  for  each  m and  n, 
for  qF  = 0 


[~— ~F  ] LJ0)(P  )K^  (q  =0)dp 

L U c+f+tA  PF  F F FF  F F 
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Equat ion  (6)  becomes  for  each  m and  n. 
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(The  higher  frequencies  are  assumed  negligible.) 

Even  in  this  simplified  problem,  a direct  so- 
lution of  the  equations  is  impracticable;  therefore 
an  iteration  procedure  must  be  devised.  It  is 
assumed  at  first  that  the  effect  of  the  after  pro- 
peller on  the  forward  propeller  (except  for  the 
thickness  effect)  is  small  and  hence  the  second 
terms  on  the  RH  of  Eqs.(l7a)and  (17b)  may  be  o- 
mi tted. 

The  iteration  procedure  is  given  below  for 
two  cases;  for  a CR  system  with  equal  number  of 
blades  (i.e.,  Nf=Na=N)  and  for  one  with  unequal 
number  (i.e.,  NF^NA) , and  neither  one  is  an  integer 
multiple  of  the  other. 


Case  if  I ; NF=NA=N 
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and  so  on,  until  values  of  loadings  are  stabilized. 
Case  til : N^Np 
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and  so  on,  until  values  of  loadings  are  stabilized. 
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It  is  to  be  noted  that,  in  contrast  to  the 
case  of  equal  number  of  blades  for  the  two  propel- 
lers of  the  CR  system,  when  the  propellers  have 
an  unequal  number  of  blades,  tyy  t Np  , the  series 
representing  the  cross-coupling  terms  due  to  the 
interaction  effects  are  more  limited.  As 
will  be  shown  later,  for  this  case  there  is  no 
unsteady  loading  Lp  on  the  forward 
propeller  at  frequency  2N-  and  no  unsteady  load- 
ing l_A  on  the  after  propeller  at  frequency  2N^. 

It  should  be  kept  in  mind  that  the  iteration  scheme 
has  been  restricted  to  the  lowest  possible  frequen- 
cies of  the  interacting  CR  system.  It  will  be 
easily  generalized  by  means  of  Eqs.(5)  and  (6)  by 
varying  parameters  ^ , £4(all  integers)  to 

other  values  than  those  already  used  (0,  ±1). 

VII.  PROPELLER  LOADING  AND  RESULTING  HYDRODYNAMIC 

FORCES  AND  MOMENTS 

Propel  I er  Load i ng 

Once  the  values  of  L^q,n^(r),  the  spanwise 
loading  components,  or  coefficients  of  the  chord- 
wise  distribution,  are  obtained,  the  spanwise  load- 
ing distribution  L^Mr)  is  determined  as  L 5 . 6 j : 


L(q)(r) 


lima*  I \ 

j iu  L^q,n,(r)9(n)sin6&deo  (19) 
o n=*l 


where  i(n)  = chordwise  modes.  Because  the  inter- 
action phenomenon  introduces  an  angle  of  attack 
even  in  the  steady-state  case  ( n)  is  taken  as  the 
complete  Bimbaum  series  which  has  the  proper  lead- 
ing edge  singularity  and  satisfies  the  Kutta  con- 
dition at  the  trailing  edge.  In  this  case  it  can 
be  shown  that 


L(q)(r) 


= l<9. 0 


(r) 


Hydrodynamic  Forces  and  Moments 

The  principal  components  of  these  forces  and 
moments  which  are  evaluated  for  each  member  of  the 
CR  system  are  listed  below  and  shown  in  Figure  2 
for  a RH  propeller  with  the  sign  convention  adopted. 

Forces:  Fx=thrust  (x-d i rect ion) 

Fy  and  Fz=hori zontal  and  vertical  com- 
ponents, respectively,  of 
the  bearing  forces 

Moments:  Qx=torque  about  the  x-axis 

(J  and  Qz=bending  moments  about  the 
^ y-  and  z-axis,  respectively 

(Subscripts  F and  A added  to  these  symbols  will 
designate  forward  and  after  propeller  cases.) 

The  elementary  forces  and  moments  can  be  de- 
termined by  resolving  the  chordwise  loadings,  act- 
ing on  ar.  elementary  radial  strip,  normal  to  the 
strip  and  taking  the  corresponding  moments  about 
any  axis  as  in  Reference  5.  However,  for  simpli- 
fication of  the  discussion  which  follows,  the 
assumption  will  be  made,  as  in  Reference  10,  that 
the  spanwise  loading  L^(r)  acts  at  the  mid- 
chord. This  is  exactly  the  case  as  far  as  thrust 
and  torque  are  concerned  and  a good  approximation 
for  the  transverse  bearing  forces  and  bending  mo- 
ments if  the  blades  are  not  wide,  which  is  the  case 
with  counterrotating  propellers. 

Thus  the  elementary  forces  acting  at  radius 
r of  an  N-bladed  propeller  will  be  given  by 


N 

1 j 

L(q»(r).' 

q(.it  + h ) 
n 

cos6  ( r) Ar 

n=»  1 

r 

N 

L 

n- I 

L(q)  ( r)  e 

iq(fjt+6n) 

s ' nSp(  r)  cos  (fit+6n)  Ar 

N 

L(q)(r)e 

iq(fit+8  ) 

n=  1 

n 

s i n6p(  r)  s i n(fit+§n)  Ar 

where  & ( r)  is  the  geometric  pitch  angle  at  ra- 
dius r and  5n=2n(n- 1 ) /N , n=l , 2 , . . . N.  The  ele- 
mentary moments  will  be  given  by 


N 

(q) 

i 

i q(fxt+5  ) 

V=i 

rL 

( r)  e 

n 

sin9p( 

r)  Ar 

N 

(q) 

i q(Qt+9  ) 

AQ  = 

y 

£ 

rL 

( r)e 

n 

cos9p( 

r)  cos (Qt+5^) Ar 

N 

-n5, 

rL 

(q) 

( r)  e 

iq(Qt+9n) 

c°s8p( 

r)  s i n(Qt+&n) Ar  (21) 

The  summation  over  all  the  blades  of  a pro- 
peller i nvol ves 

1)  for  thrust  and  torque 

N i q6  rN  when  q=nN,  n=0,l,2  ... 

e n = { (22) 

n=!  '•0  for  all  other  q 

2)  for  transverse  forces  and  bending  moments 


N i(q±!)§ 


n when  q=nN+l , n=0, I , 2 ... 
'■0  for  a 1 1 other  q 


It  is  thus  evident  that  in  the  steady-state 
case  (q=0)  thrust  and  torque  from  each  propeller 
will  be  present  (see  Eq.22),  whereas,  since  the 
condition  under  consideration  is  that  of  uniform 
inflow  into  the  forward  propeller,  there  will  be  no 
transverse  forces  and  bending  moments  (see  Eq.23). 
This  will  be  so  whether  the  propellers  are  of  equal 
or  unequal  number  of  blades. 

As  has  been  shown  the  interaction  phenomenon 
induces  unsteady  loadings  on  both  units  of  the  CR 
system.  In  the  case  of  equal  blade  number, those  on 
the  forward  propeller  are  Lp  '(rp),  on  the  after 

propeller,  L^^(r  ).  As  seen  from  Eqs.(22)  and 
(23)  both  propellers  generate  thrust  and  torque  at 

frequencies  q = 2jtN,i=l,2 which  corresponds 

to  blade-blade  crossing  frequency  for  such  propel- 
lers, and  no  unsteady  transverse  bearing  forces 
and  moments,  s i nee  no  combination  of  integers  £ and 
n can  satisfy  the  relation  (2£-n)N=JI. 

In  the  case  of  unequal  blade  number,  N^Np, 
the  unsteady  loadings  on  the  forward  propeller  are 
at  frequencies  qp^^N^  and  2£aNp  and  on  the 
after  propeller  at  q^=2^3N(:  and  2£4N/\  . The 

criterion  for  thrust  and  torque,  Eq.(22),  yields 

qp  = " "iNF 

qF  = 2£aNp  = n-Np 

qA  = 2£aNF  = "A 

\ " 2^NA  = nA  (22a) 

and  the  criterion  for  transverse  forces  and  bending 
moments,  Eq.(23),  yields 
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„r  = ?-'^A  - njNf.  + 1 
a - / i>snf  'I 

\ " 2 ,,Mr  " n'NA  Tl 

nA  ” ?C‘MA  y "A  (23a) 


The  condl  t ions  of  ( 2 2 .t  ) for  thrust  and  torque  arc 
•i  1 ways  satisfied  for  frequency  q’qF=qA  by  choos- 
ing j\  * iiitlr,  fa’  mllA,  m-1,2,  ...  , so  that 
q - ? , i.e.,  the  so-called  blade-blade  cross- 

inq  f ri  qtrrcy  and  multiples  thereof.  Cqua  1 1 on  (22a) 
r in  iho  he  satisfied  at  lower  frequencies  (I)  if 
Nv"kl.'.  , '<  heinq  an  inteqer,  in  which  Case  choosing 
fjk  e 2;,  » I,  " (,  k yields  q " 7mkNe , and  (2)  If 
mil  |.'|  are  both  even  numbers,  In  wliich  case 
tl  «■  m'-.f.y.  In  tie  latter  two  cases  the  conditions 
nf  (2)il  far  transverse  fo rats  and  bending  moments 
are  obviously  not  satisfied. 

When  MA^kNF  and  UA  and  Up  are  not  both 
even  nin'm1.,  the  C it  system  generates  thrust  and 
torque  only  at  q " 2niV«Np  (blade-blade  crossing 
I ri  queer,  i ns)  . 

The  conditions  of  (23a)  arc  satisfied  for 


q - qp  +'  » 2?jHa  + ' “ njHp 
')  ” 1A  ±l  " 2/jl'F  " n3Nft 


since  2 - n^N  +n  1^ 

•jiH'ocl'is  *v.  v/l»  i <■  1 1 side  forces 
of  I he  CR  systeni  occur  are 


Is  possible.  The  f re- 
am! bending  moments 


q «■  2 fj  N Tl  •>  2 f3nr  ±1 

from  which  ? ^ N«  *s  2 /.3Nr  t2  . 

An  • c,sy  way  of  determining  the  frequencies 
for  al'eriwt:nq  forces  and  moments  Is  to  write  out 
the  s».  cyuence  of  numbers  which  are  integer  multiples 
of  twice  the  blade  number.  for  example,  for  3 sod 
5 blades 


<°*V 


-Nr  e 
I o 


I q(Jt 


-N.  r e 
A o 


i qflt 


jL^fr^sinOptrp)  rpdrF 
o 

lLiC')<rA)sinVrA)V,rA 

o 

(2l,b) 


where  q’O  for  steady  state,  and  the  lowest  frequency 
of  the  alternating  thrust  and  torque  is 

q <*  2kNp  when  N^kNp,  k"l,2,  ... 

q « N^Np.  when  tl^Np  end  both  are  even  num- 

q «=  2NAr,F  f°r  al'  other  Nft/Np 


2)  for  transverse  bearing  forces  and  bending  mo- 
ments. wlien  N.t*kU  and  N.  and  N_  are  not 
. i A F A F 

both  even  numbers. 


'Vf-T 


F 

sr  r c 
? o 


i (q-+1)ttt  I (qF) 

Jlf  F(rF)sinep(rF)drF 


(F7)F  = ?i(Fy)F 


Na  l(q*t)nt  I (qA) 

I^a’TV  ILA  (rA)5inUp(rA)drA 

' O 


<Fr>A-*i(,VA 


(25a) 


tl,  , Kq-^OOt  I,  (qj 

(VF’Tro  C JLF  (rF)c05yrF>rF<lrF 

' o 


3:  6 12  If  24  30  ... 
5:  10  20  30  ... 


The  frequencies  lor  side  forces  will  be  tile  mean  of 
any  pair  ol  numbers  In  the  two  sequences  which  differ 
by  2,  in  this  case  !!,  19  ...  . The  frequencies 
for  thrust  aid  torque  will  hr  the  mean  ol  any  pair 
ol  numbers  which  are  alike,  in  this  case  30, 6o...  . 
for  6 end  *1  b lades,  the  two  lowest  frequencies  will 
be  2b  and  hfi  for  thrust  and  torque.  (There  is  no 
side  force  in  this  case.) 

Another  derivation  of  the  frequencies  of  the 
alternating  lories  developed  by  interactions  be- 
tween a nair  of  counterrot  at Inq  propellers  in  a 
uniform  :nflow  is  O'ven  by  Reference  II,  with  the 
sane  results. 

On  the  basis  of  the  preccdinq  discussion,  t he 
total  forces  and  moments  ere  obtained  from  (20)  and 
(21)  as; 

I)  for  thrust  and  torque 


^Fx3f" 

<F*V 


N_.r  r 
V c 


N.  r e 
A o 


IqOt 


i qUt 


f LF'’)(rF)cos0p(rF)drF 
o 

J’  LAq)(rA)c05VrA>drA 


(2ba) 


«jz)F*  ;.«yF 


, A a 

(Va"—  ro  c 


i(qAll)Qt  I (q  ) 

JLA  (rA)cOS0P(rA)rA<JrA 


(QZ)A“  *'<Va 


(25b) 


where  the  upper  signs  (at  q"qp- 1 trqA't  I ) are  used 
when  Na  > NF  and  tlie  lower  siqns  (at  q°qF+  I =q.  - 1 ) 
are  used  when  NA  Hp,  and  qF=2^NA  and  qA’2i3llF 
where  Vj  and  i3  must  satisfy  the  condition 
2t1NA-2t3NF*2  . 

It  is  to  he  noted  l ha t r0“rF(j»  forward  pro- 
peller radius  and  rA  and  rp  are  fractions  of 
rFo,  and  that  finally  the  real  parts  of  the  forces 
and  moments  are  to  he  taken. 

Blade  Pending  Moment. 

Following  Reference  6,  the  blade  bending  mo- 
ment about  tie  face  pitch  line  at  any  radius  r. 
of  a blade  of  the  forward  or  after  propeller  is1 
calculated  from  the  spanwise  loading  L'c*'(r)  at  any 
shaft  frequency  q as 

m(9)  =r3  c'qOt  j L(<,)(f.)c05[  0p(r)  - 0p(  rj)  ] ( r-r. ) dr  (26) 


The  instantaneous  blade  bending  moment  distribution 
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! 

[ 


a 


as  t Hr  blades  sw'nq  around  the  shaft  Is 

M - Re^Hh(q)e,q0t  (2/) 

b (|  b 

Here  q-Hj.  for  the  forward  propeller  and  q=q^  for 
tie  af  ter  propel  ler. 

VIII.  NUMERICAL  RESULTS 

On  the  ha*.  Is  of  the  theoretical  procedures 
outlined  i rs  the  preceding  sections,  a numerical  ap- 
proach has  heeii  established  and  adapted  to  the 
COC-f’f’OO  or  7 V)0  high-speed  digital  computer.  The 
program  furnishes  in  the  case  of  uniform  inflow  to 
the  forward  propeller,  a)  the  steady  and  time-de- 
pendent blade  loadings,  h)  the  corresponding  hydro- 
dynamic  forces  and  moments,  and  c)  the  Made  bend- 
ing moment  al-orr*  the  face-pitch  1 ine  at  any  radius, 

T lie  p>  press  Ions  for  tile  kernel  functions 
qiven  by  Eg, .(7-9),  those  for  the  normal  velocities 
on  the  left -hand  sides  of  the  integral  equations 
given  by  Eqs.  (11-16),  and  the.  pair  of  integral  e- 
qualions  (l/a,h  and  I8n,b)  , constitute  the  desired 
working  forms.  The  computer  program  prepares  all 
the  necessary  inlormation  for  the  execution  of  the 
suggested  iteration  procedure.  It  is  a lengthy 
ptogram,  tin:  duration  depending  on  tire  number  of 
blades  of  each  component  of  the  CR  system,  on  the 
nurber  of  selected  chordwise  inodes,  and  the  number 
of  frequencies  for  which  information  will  be  dc- 
s I red. 

To  establish  the  accuracy  and  usefulness  of 
the  computal  ion.il  procedure,  a correlation  with 
(Misting  experimental  results  must  be  made.  Calcu- 
lations have  been  performed  for  two  CR  configura- 
tions for  whici  data  are  available  from  tests  at 
David  Taylor  NS R DC  l 12, 13].  The  data  are  for  the 
two  lowest  frequencies,  and  fj-1,  and  the  com- 

putations have  been  limited  lo  these  frequencies. 

Hu  iterative  procedure  starts  after  all  the 
required  information  has  been  computed  and  stored 
properly.  The  necessary  kernel  functions  and  in- 
verse kernel  functions  are  calculated  for  both  a 
b-O-b  ant  .1  b-0-5  CR  system  for  values  of  m and  q 
as  indicated  in  the  table  below. 

TABLE  I 

b-O-b  CR  System  b-0-5  CR  System 


indicates  that  the  number  of  iterations  need  not  be 
greater  than  10,  although  the  execution  time  of  an 
additional  iteration  is  minimal. 

The  b- 0-b  CR  system  is  composed  of  the  Oovid 
Taylor  IJSRDC  propeller  3686  forward  and  propeller 
3687-A  aft,  and  the  **-0-5  set  of  propeller  3686 
forward  and  propeller  3869  aft  L 1 3 J - Propeller 
characteristics  an-!  flow  conditions  are  qiven  in 
the  fol lowi ng  table. 

TABLE  2 


Propeller  Propeller  Propeller 

3686  3687-A  38**9 


Number  of  Blades.N 

b 

b 

5 

EAR 

0.303 

0.322 

0.379 

P/D  at  0.7R 

1.291 

1.320 

1.287 

Diameter,  D,  in 

12.017 

11.776 

1 1 . 785 

Rotation 

L.H." 

it 

R.H. 

R.H." 

n,  rps 

12 

12 

12 

Speed,  ft/sec 

13.22 

13.22 

13.22 

Advance  rate,  J 

l.l** 

l.l** 

1 . 1 

>’•'{_ . H . rotation  Is 
R.H.  rotation  is 

ccw  looking 
cw  looking 

forward 

forward 

**J  based  on  diameter  of  forward  propeller  36cb 


The  theoretical  predictions  with  and  without 
thickness  effects  and  the  experimental  results  are 
summarized  in  the  following  table. 

TABLE  3 

CORRELATION  0E  PREDICTED  AND  MEASURED  VALUES 
b-O-b  CR  SYSTEM 

Forward  Propel ler  3686 
After  Propeller  3687-A 

Axial  Spacing  " 0.28  forward  propeller  radius 
J-I.l 


: 


I 


j 


I 


i 


Inverse  Inverse 

Matrix  Matrix  m q Matrix  Matrix  nr  q 


Q = 0 , Steady  State 


EXPERIMENT 


AF 

kaf 


AF 

kaf 


'fa 

Va 


'FA 

Va 


FF 


FF 


AA 


AA 


0 0 

*l  0 


b 8 
8 8 
-8 

0 0 
b 0 


b 8 
8 8 
8 


AF 


AF 


'FA 


'FA 


0 

0 

W i thout. 

With 

Thickness 

Thickness 

Vf 

0 

Forward 

<VF 

0.126* 

0.1639* 

0.125 

5 

10 

Forward 

(kVf 

0.0295" 

0.0373* 

0.0315 

After 

<Va 

0. 1b 5* 

0.1335* 

0. 150 

kff 

-10 

After 

(Va 

0.0330* 

0.0308* 

0.0315 

0 

0 • 

0 “ 2N  ” 

8 

K . 

0 

Forward 

(Vf 

0.0b85 

0.07b7 

0.0285 

AA 

0.0058 

Forward 

(KCV 

0.0101 

0.015b 

b 

8 

0.152b 

After 

<Va 

0.1156 

0.0095 

< 

< 

8 

After 

0 . 02**2 

0.0312 

0.0022 

A workable  program  has  been  established  which 


Including  friction 
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CORRELATION  OF  PREDICTED  AND  MEASURED  VALUES 
<1-0-5  CR  SYSTEM 

Forward  Propeller  3686 
After  Propeller  3869 

, , Axla|  Spacl  ng=0. 28  forward  propeller  radius 
J-  1.1 

Q a 0,  Steady  State 

THEORY EXPERIMENT 


W i thout 

With 

Thickness 

Th i ckness 

Forward 

<Vf 

0.087* 

0 . 1 060* 

0.130 

Forward 

(kVf 

0.0207* 

0.0266" 

0.030 

After 

<*T>A 

0.161* 

0.1093* 

0.130 

After 

<*Q>A 

0.0361* 

0.0255* 

0.028 

Ip  = ?HA 

-1=9 

„ 

forward  (1^)^ 

0.0079 

0.0239 

0.0075 

Forward 

<*FV>F 

0.0079 

0.0239 

0.0076 

F orward 

^QflV 

0.0062 

0.0053 

0.0060 

Forward 

(i?'QvV 

0.0062 

0.0053 

0.0061 

*a  - 2Nr 

.+1=9 

After 

^Fli^A 

0.0231 

0.0275 

0.0057 

After 

<Rrv}A 

0.0231 

0.0275 

0.0066 

After 

0.0159 

0.0181 

0.0023 

After 

^QV^A 

0.0159 

0 . 0 1 8 1 

0.0023 

if 

Including  friction 


IX.  CONCLUSIONS 

Linearized  unsteady  lifting-surface  theory 
has  been  applied  in  the  study  of  two  Interacting 
propellers  of  a counterrotating  system,  when  both 
, units  operate  in  a spatially  non-uniform  inflow. 

A mathematical  model  Is  introduced  taking  Into 
account  the  exact  geometry  of  the  propulsive  sys- 
tem as  well  as  the  three-dimensional  spatially 
varying  inflow.  The  propeller  blades  are  con- 
sidered to  be  of  finite  thickness  and  lying  on  a 
helicoidal  surface  of  varying  pitch.  The  blades 
have  arbitrary  planform,  camber  and  sweep  angle. 

The  computational  procedure,  however,  has 
ecn  developed  and  adapted  to  the  COC  6600  and 
*7600  high-speed  digital  computer,  for  the  case 
where  both  units  operate  with  the  same  RPM  In  a 
uniform  inflow  field. 

'The  uniformity  of  the  inflow  field  provides 
for  a be t ter ' unders tand I ng  of  the  mechanism  of 
Interaction  of  the  CR  system,  since  t lie  presence 
of  wake  harmonics  would  have  such  a dominant  ef- 
fect as  to  mask  the  interaction  phenomenon. 

The  study  provides  information  about  the 
Steady  and  unsteady  blade  loading  distributions 
and  the  corresponding  hydrodynamic  forces  and  mo- 
ments on  both  components  of  the  propulsive  device. 

Rules  have  been  established  for  the  presence 
or  absence  of  the  steady  and  unsteady  hydrodynamic 
forces  and  moment  when  the  CR  system  is  made,  up 
of  equal  and  unequal  numbers  of  blades.  In  fact, 
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when  the  propellers  of  the  CR  system  have  equal 
number  of  blades,  only  the  steady  and  unsteady 
thrust  and  torque  will  be  generated  on  each  pro- 
peller of  the.  CR  system,  at  zero  and  blade-blade 
crossing  frequencies  or  multiples  thereof  (q=2/N, 

2=1, 2, 3 ...  and  N~common  number  of  blades).  When 
the  CR  system  is  made  up  of  propellers  with  un- 
equal number  of  blades,  i.c.,  Nr/N^,  there  is  an 
easy  way  of  determining  the  frequencies  of  al- 
ternating (unsteady)  forces  and  moments. 

The  numerical  work  has  been  limited.  The 
calculations  without  thickness  effect  have  shown 
a reverse  trend  to  that  of  experiment  for  the  un- 
steady forces  and  moments  generated  on  the  after 
propeller.  The  numerical  calculations  have  shown 
much  higher  values  for  the  forces  and  moments  ex- 
erted on  the  after  propeller  than  those  on  the 
forward  propeller.  It  is  difficult  to  explain 
the  discrepancy  until  more  detailed  and  exhaustive 
numerical  calculations  are  performed  (including 
thickness  effect).  Furthermore,  additional  ex- 
perimental work  should  also  be  conducted. 

Some  remarks  arc  in  order  on  the  discrep- 
ancies between  theory  and  experiment. 

1)  A few  years  ago  an  attempt  was  made  to  cal- 
culate unsteady  thrust  and  torque  of  a 
counterrotating  propeller  system  by  a ■•trip- 
wise  method  following  the  Kemp  and  Sears 
approach  for  the  mutual  interference  of  blades 
in  cascade  due  to  relative  motion  of  succes- 
sive blade  rows . Results  of  this  study  have 
shown  that  the  unsteady  thrust  of  the  alter 
propeller  is  much  higher  than  that  of  the 
forward  propeller,  when  the  bound  vorticity 
and  the  vortex  wake  shed  by  the  forward  pro- 
peller are  taken  into  account,  either  with 
flat  plate  inodes  or  with  an  elliptic  distri- 
bution. 

2)  Furthermore,  observation  of  Dr.  Wereldsma's 
experimental  results  (although  these  are  for 
non-uniform  inflow  conditions),  especially  of 
the  time  records  of  the  horizontal  and  verti- 
cal blade  bending  moments,  show  that  super- 
imposed on  the  fluctuations  due  to  the  wake 
at  blade  frequency,  are  the  fluctuations  due 
to  blade  interference  at  blade-blade  crossing 
frequency.  It  can  be  seen  that  these  fluctu- 
ations due  to  Interference  arc  much  larger 
for  the  after  propeller  than  for  the  forward. 

The  question  of  the  discrepancies  thus  cannot  be 
considered  a closed  subject,  but  requires  further 
Investigation  both  theoretical  and  experimental. 
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:ig.  1 : Counterrotating  propeller  arrangement  - 

angular  coordinates. 
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Fig.  2:  Resolution  of  forces  and  moments. 
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APPENDIX  A 

Effect  of  Blade  Thickness  of  Each  Propeller 
on  the  Velocity  Field  of  the  Other 


I.  The  thickness  distribution  of  a blade  section  is  represented  by  a 
source-sink  distribution  assumed  to  be  smeared  over  a projection  of  the 
section  in  the  propeller  plane.  The  velocity  potential  due  to  blade  thick- 
ness of  the  after  propeller  at  a point  (xp,  Tp,cPp)  on  the  forward  propeller 
is  then  given  by 


(*p(xp,rp,cpF;t); 


N.  0UA 

A bA 


TTn  * C I 


n=i  -ew  pa 


A A AO 


/ ffi  n \ 2U&f(  §A,  PA,  ®Aq) 

where  M(  PA, 0AO)  = — source  strength  density  determined 

"A 

in  accordance  with  the  "thin  body"  approach, 


^ ^A*  ®A(P  = th'c'<ness  distribution  over  one  side  of  the  blade 

section  at  radial  distance  PA  in  the  propeller  plane 

RAF  - (<V5i>2  + 1 + "A  ■ SrFPA  “^“ao^FO-20'*^3}'72 

XF  " W"  “ (oF  - 8bF  costfa  /a 

V - V8  + * * <°A  - 9bA  “W'  + 6 
0A  = (2ti/Na)  (n-l),  n = 1,2,...NA 


Since 


SbA  sin°a  ^ 


Eq.  (A-l)  can  be  reduced  to 


, , u Na  f"  . sf<<Ve»> 

<sf>  - - A * J J — ST — -nr-*-  »•. 


n=l  o P« 


A-l 
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The  thickness  distribution  f(  PA>  6a)  will  be  approximated  by  a 
lenticular  section,  i.e., 

t(  p 


f(pA’0a)~  “5^  Sin2ea 


— ( p ) • pA  6 sin2 6 
c A ' A bA  or 


(A-3) 


where  T is  maximum  thickness  in  the  projected  plane,  t J c is  ratio  of 
maximum  thickness  to  chord  of  the  expanded  section,  and  is  projected 

semichord.  In  addition,  the  reciprocal  of  the  Descartes  distance  R^.  is 
expanded  as 


J_  = I z e'me  J (IK)  e A dk 

l\a  TT  m „ m 

AF  m=— 00  -°° 


(A-M 


where  0 . + <PF0  - SDt  + eAn 


and  ( I K)  = 
m 


lm(lklPA)  Km(lklrF)  for  PA  < rF 
'm('kl  rF)  Km(lklPA)  for  > rp 


N 


A imB, 


With  I e 
n=  I 


An 


Na  if  m = X Na,  X = 0,  ±1,  -2, . . . 


0 otherwise 


and  If  ~ 2 -r  ( Pa)  • PAebA  Sin0o  COS  eo 

or 

the  velocity  potential  ( A— 2 ) becomes 

UN.  “ -O  n<.  im(9rn+CTA) 

(9  ) w-  — Z e-'2mat  e F0  A 

F tA  tt2  m=-® 
m=XN. 


11  - im6,  .cos 0^  t 

I J e “ s i n cos  0bA  -2(  pA ) J)  +a2  p2A 

O Pa 


J 


'(vP1 


dk  dPA  decr 


(A-5) 


A -2 
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The  velocity  normal  to  the  forward  propeller  blade  due  to  the 
thickness  of  the  after  propeller,  which  is  an  additional  perturbation 
on  the  L.H.  of  the  first  integral  equation, is  derived  as 


(Wp)tA  ' " 3np  ^tA  “ ’ Vl+a2r2F  ^ dxF  " r2f  ^ F't A 


Therefore, 


(^) 


irFNA 


tA  T^Vl+a2^  m=-°° 
m=  IN. 


v -i2mQt 
L.  e e 


im(aA+oF)  -im6bf.  coscp^ 


-im0,  . cos  0 
p bA  a 


Sin0a  coS6(<  CA6bA  -S.  (pA)/TW^ 


I-”,...,  J(W8e)k/>  »' 

J (|K)me  Kak  ~ ~r)  e dkdpAdec 


(A-6) 


The  0 - integral  involves 


71  i(  k/a-m)  0,  .cos  0 
r v ' ' bA  a 


sin©  cos0  d 0 
a a a 


? i(u/a)0bAcosea 

= J e sin0^  cos©^  d©&  ( letting  u=k-am) 


F(u, P.) 


(A-7) 


r 5in<uWa)  - (uVa)  c°5<uVa>i 

where  F(u,  p.)  = 1 


After  substituting  (A-7)  and  u=(k-am)  in  (A-6)  and  applying  the 
generalized  lift  operator  (see  Eq.  7a),  the  velocity  becomes,  for  each 
lift  operator  mode  m , 
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2a2rJ, 


® -i2mQt  im(2a  -ae) 
£ e e 


tA  r^Vl+a2^  m=- 00 
m=*£N^ 

f PA  t0<PA) 

■ J T„  ~ 

p bA 

a 

09  i(a_-a  -ae)u/a 

* J I m(  1 u+am I p . ) K ( I u+aml rc)e  (au+a2m-  ~) 

* fn  n rn  r p 

rF 


I m ((-u-2am) 0bF/a)  F(u, PA)  du  dpft 


(A-8) 


for  PA  < r^,  otherwise  pA  and  are  interchanged  in  the  modified  Bessel 
functions.  On  changing  the  doubly  infinite  u-integral  to  an  integral 
from  0 to  +“  , (A-8)  can  be  written  as 


V u 4a 


g,VA  ~ TSmtlt  -'"(2Vae) 

- L & e 


tA  T^yi+a^r2  m=-°° 
y m=jKNA 


OB 

* l T"“T  (PA^  ^ 1+a2p2A  I f<u'pa> 
PA  bA  o 


r i(c  -a  -ae)u/a 

• I (lu-am|p.)  K (lu-amlr_)  e 

• (au-a2m+ l^m4  ( -u+2am)  0,  /a) 
rp 

-i(a  -a  -ae)u/a 

+ I (l u+aml pA)  K ( I u+aml r)  e 
m Am  r 

. (au+a2m-  (u+2am) 0bF/a)  L dudPA  (A-9) 

rF 

where  the  sign  of  m is  also  changed,  which  is  permissible  since  the  m series 
is  doubly  infinite.  The  integrand  of  (A-9)  is  zero  when  u=0  since 
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F(0, Pfl)  = l im 

u-o 


sin(u9  /a)  - (ue  /a)  cos(u0  /a)' 


In  the  steady-state  condition,  n*=4=0,  the  velocity  on  the  forward 
propeller  due  to  the  blade  thickness  of  the  after  propeller. 


^°’a)(rF)j 


can  be  shown  to  be  as  given  in  Eq.  (13)  of  the  text.  The 


unsteady  velocity  components  are  at  frequencies  2mfi  = 24n  f},  with  4 = -1,-2,., 
and  are  given  by 


(o  (m,m)\ 

■h-L  *' 


/ ( I AIN  ,m)\ 


/.  (-U1N  ,m)\ 

|WF  J -i2|i|NA0t 


8 (|1|VS)\  (0 

+ conj  I 


12 l 4IN  Ot 


/ (UIN  ,S)\ 

/W_  J i2|4|N.0t 

■ s{^r-L  * (fl-,0) 

(a  (,*"VS)\ 

Herel— — I is  the  coefficient  of  exp( i2mftt)  in  (A-9) 

\ U /tA 

with  m = 141Na  and  is  half  the  velocity  at  2|4|Na  times  the  shaft  frequency. 
This  velocity  is  designated  as 


. (24N  ,m) 

WF  <rF) 


and  is  given  by  Eq.  (14)  of  the  text. 
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II.  The  velocity  potential  due  to  the  blade  thickness  of  the  forward 
propeller  at  a point  ( xA,  r/\>  on  the  after  propeller  is  given  by 
(cf.  (A-l ) and  (A-2)) 

TT 


(VVa'V1'. 


U / "r  WV 

2V  \ J J — Se 


V,A,TA,LVke  “ " 2-n  “ J J B1T~  R_. 

tF  n=l  o Pp  a FA 


2 2 2 


dPF  d6o 


(A-l 1) 


where  RpA  = + rA  + Pp  - 2rft PpCOsL- VW^+Vn J; 


1/2 


XA  = V3  + 6 = (CTA  - ®bA  cosV/3  + 6 
V = Wa  = (°F  " 6bF  cosV8 


and  again  approximating  the  thickness  distribution  by  a lenticular  section 


Bf 

a 


PF  ®bF  sinVos0Q 


The  velocity  normal  to  the  blades  of  the  after  propeller  due  to  this  velocity 
potential  is 


<V 


tF 


( A— 12 ) 


If  the  procedure  of  the  preceding  section  is  followed,  it  can  be 
easily  shown  that  the  velocity  normal  to  the  blades  of  the  after  propeller 
induced  by  the  forward  propeller  thickness  is  in  the  steady-state  given  by 
Eq.  (15)  of  the  text  and  in  the  unsteady  case  at  frequencies  equal  to 
2 1 X | Np.  times  the  shaft  frequency  given  by  Eq.  (16). 


